Actively contracting bundles of polar filaments 
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We introduce a phenomenological model to study the properties of bundles of polar filaments 
which interact via active elements. The stability of the homogeneous state, the attractors of the 
dynamics in the unstable regime and the tensile stress generated in the bundle are discussed. We find 
that the interaction of parallel filaments can induce unstable behavior and is responsible for active 
contraction and tension in the bundle. Interaction between antiparallel filaments leads to filament 
sorting. Our model could apply to simple contractile structures in cells such as stress fibers. 
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Living cells have remarkable mechanical properties. In 
addition to a passive response to mechanical stresses, eu- 
caryotic cells are able to actively change their shapes, to 
generate motion and forces Q and to react to externally 
imposed mechanical conditions H. The cytoskeleton, 
which is a complex network of elastic protein filaments 
such as actin filaments and microtubules, plays a key 
role in these processes. Filaments are rod-like structures 
that interact with a large number of specific proteins || . 
Examples are cross-linking proteins which induce forma- 
tion of a gel-like filament network, and bundling proteins 
which lead to filaments aligned in parallel. Furthermore, 
motor proteins of the cytoskeleton are able to use the 
chemical energy of the hydrolysis of Adenosinctriphos- 
phate (ATP) to generate forces and motion along fila- 
ments Myosins for example interact with actin 
filaments while kinesins and dyneins move along micro- 
tubules. The direction of motion of motor molecules is 
determined by the polar structure of the filaments which 
have two different ends, one denoted "plus", the other 
"minus" : a given type of motors always moves towards 
the same end . Small aggregates of motors which con- 
tain two or more active domains can bind at the same 
time to two filaments and exert relative forces and motion 
between them. In contrast to passive filament solutions 
whose rheological properties have been studied in recent 
years Q , filament systems which interact with molecular 
motors represent intrinsically active materials which ex- 
hibit rich types of behavior that have recently attracted 
much interest. In particular, the formation of asters and 
spiral defects, as well as the shortening of filament bun- 
dles has been observed in vitro and in vivo §0. The 
self-organization of motor-filament systems has also been 
addressed theoretically [pl]dj3[ . 

Bundles of actin filaments interacting with myosin mo- 
tors are prominent cytoskeletal structures that are in- 
volved in many active phenomena in the cell For 
example, within the sarcomeres of skeletal muscle fibers, 
they are responsible for muscle contraction; as stress 
fibers they produce forces in cells; and as contractile rings 
they are important during the final step of cell division. 
The contraction of sarcomeres follows from a particular 
arrangement of myosin and actin filaments. In contrast, 
stress fibers and related structures lack an obvious spatial 



organization of their components. The existence of these 
simple contractile structures raises interesting questions: 
What are the minimal requirements for a bundle of fila- 
ments and active elements to contract? Do stable steady 
states of the bundle exist which generate tension? 

In this Letter, we develop a simple phenomenologi- 
cal model for bundles of polar filaments containing ac- 
tive elements, that allows us to answer such questions 
and to estimate the mechanical tension generated in the 
bundle. Furthermore, we perform computer simulations 
which demonstrate that the general behaviors described 
by this model are indeed obtained for simple realizations 
of the motor-filament interactions. Finally, we relate our 
findings to biological systems. 

Consider a linear bundle of aligned polar filaments of 
equal length £, distributed along the rc-axis. The number 
densities of filaments with their plus-ends pointing to the 
right and to the left and with their center located at 
position x are denoted c + (x) and c~(x), resp. These 
densities satisfy the conservation laws 



d t c + = Dd 2 x c^ 



d t c~ = Ddtc- - d x J- + - d x J~ 
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FIG. 1. Schematic representation of motor-filament inter- 
actions. Motors are assumed to move towards the plus-end 
of filaments, arrows indicate the direction of filament motion. 
Vertical lines indicate the centers of gravity of the filament 
pair, (a) Antiparallel filaments slide in opposite directions, 
(b) Relative motion of parallel filaments occurs if a motor 
binds to a plus-end. Parallel filaments thus have the ten- 
dency to align their plus-ends, (c) Tension profile a along a 
filament driven by a point-force. 
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Here, D is an effective coefficient for filament diffusion 
along the ir-axis. The currents J H and J h are ac- 
tive filament currents that result from interactions me- 
diated by motors between pairs of antiparallel filaments. 
Currents resulting from interactions between parallel fil- 
aments are denoted J ++ and J . Here, we have as- 
sumed for simplicity that two-filament interactions dom- 
inate, which corresponds to a sufficently low density of 
active elements. Eq. (f) thus corresponds to a coarse- 
grained description where details of the dynamics of the 
motors have been eliminated. The direction of filament 
motion induced by a motor only depends on the fila- 
ments' relative orientation, see Fig. 1. Therefore, we can 
use general symmetry arguments to write expressions for 
the active currents without referring to a specific interac- 
tion mechanism. To this end we express the currents in 
terms of the filament densities c + and c~ . For example, 
we write J H (x) = f_fd4 j H (x,£), where j H (x,£) is 
the average current of "minus" -filaments located at x in- 
duced by the interaction with "plus" -filaments located at 
x + £. The integral arises since all filaments with centers 
located in the interval [x — £,x + £] can interact with a 
filament located at x. Then, assuming that the prob- 
ability for an active interaction between two filaments 
only depends on the probability of their encounter, we 
write j H = v + ~ (£)c + (x)c~ (x + £). Here, w H (£) 

is the effective relative velocity between two antiparal- 
lel filaments a distance £ apart induced by many indi- 
vidual events of motor activity. For parallel filaments 
we get analogously, e.g., j ++ = v ++ (£)c + (x)c + (x + £). 
Now, symmetry arguments restrict the possible form of 
the functions w ±:t and w ±=F : In the absence of exter- 
nal forces momentum conservation requires that the cen- 
ter of gravity must remain fixed when two filaments 
are displaced. This leads to i> ±=t (£) = — v ±:t (— £) and 
(£) = — v h (— £), see Fig. 1. Furthermore, invari- 
ance of the system with respect to inversions of space 
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FIG. 2. The critical value a c = a c £ 2 c/D as a function of 
/3 — /3£cL/D for different values of L = L/l. Here, a and /3 
characterize the interaction strength of parallel and antipar- 
allel filaments, respectively, L is the system size and £ the fil- 
ament length. Solid lines correspond to 8c — c/10 and dashed 
lines to Sc = c/2, where c is the filament density and Sc/c the 
relative difference of the densities of plus and minus filaments. 
The homogeneous state is unstable for a > ct c . 



leads to the condition = — v (— £). Respect- 

ing these criteria, we choose for simplicity v '"(£) = 
— u H (£) = (3 and = «~~(£) = a sign(£). Here, 

a and (3 are constants and sign(£) = ±1, depending on 
the sign of £. We have checked that other choices do not 
alter the general properties of our model. The currents 
thus read 



(2) 



J ±± (x)=a / d£ [c k (x + ^) -c ± (x-^)]c ± (x) 
Jo 

r e 

J±t ( x ) = T /3 d£ c T (x + £)c ± (x) . 
J-t 



(3) 



Equations (Q) and @ describe the dynamics of our 
model. We first analyze the stability of the homo- 
geneous state with constant c ± (x) = cj, which is a 
fixed-point of the dynamics for all values of the pa- 
rameters. Using a Fourier expansion of the densities 
c ± (x) = Cq +J2k c^e lkx , the filament dynamics to linear 
order in the amplitudes is given by 

Af4 Wa++ a+-\(4 

dt\cj-)-\A-+ A-){c k 

where the elements of the matrix A(fc) are 

A ±± (fc) = -Dk 2 - 2a(cos(fc£) - l)cfi ± 2ij3klc^ 
A ±T (fc) = ±2if3sm{k£)c Q t . (4) 

With A(fc) denoting the larger of the real parts of the two 
eigenvalues of the matrix A(fc), the homogeneous state is 
stable if A(fc) < for all k. When an instability occurs, 
a band of unstable modes appears, which extends from 
k = to some positive k. For a system of size L with 
periodic boundary conditions (e.g., a contractile ring), 
the stability of the homogeneous state is determined by 
the sign of A(fc m ; n ), with fc m j n = 2ir/L. We find that the 
homogeneous state becomes linearly unstable as soon as 
a > a r with the critical value 



D_ (£3£cL_ Sc L 
Pc * \ D ' c' £ 



(5) 



Here c = + c$ , Sc = Cq — Cq , and f(u,v,w) is a 
dimensionless scaling function. Corresponding stability 
diagrams are displayed in Fig. 2. For all values of the 
parameters we find a c > 0. Therefore, unstable behavior 
is induced by an interaction of parallel filaments with 
a > 0. Consistent with this finding a c decreases with 
increasing Sc, i.e., with an increasing fraction of filaments 
pointing in the same direction the system becomes less 
stable. The critical value a c increases with increasing D, 
indicating that filament diffusion has a stabilizing effect. 
In many limiting cases simple analytic expressions for a c 
can be obtained. In particular, we find in the limit of 
large L, f = 2/(1 + \Sc\/c) if f3 = or Sc = and / = I 
otherwise. In the homogeneous state, the total filament 
current J = J ++ + J H + J h + J~~ vanishes, while 
currents of plus and minus filaments exist in opposite 
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FIG. 3. Variance of stable stationary filament distributions 
as a function of a = a£ 2 c/D obtained by numerical solution 
of Eqs. (1) and (2) for L = 10£ and a system containing 
plus-filaments only (c~ = 0). The homogeneous state is sta- 
ble for a < d c . For a > fij a new attractor corresponding 
to a shortened bundle appears. A third attractor is also in- 
dicated. The insets show for each attractor typical filament 
distributions over one period in arbitrary units. 

directions. As soon as the homogeneous state becomes 
unstable, a transient current J occurs which redistributes 
filaments until an inhomogeneous stationary state is at- 
tained. This final state consists of one or several short- 
ened filament bundles separated by empty regions. Fig. 3 
displays the variance < x 2 >= j dx x 2 (c + + c~) of 
stationary filament distributions as a function of a, 
where we have assumed without loss of generality that 
J dx x(c + + c~) =0. For small a only the homogeneous 
state exists as attractor of the system. A second attractor 
appears for a = ay < a c , which corresponds to a short- 
ened filament bundle. A third type of attractor obtained 
numerically consists of two shortened filament bundles. 
However, we cannot distinguish numerically whether this 
state is a real attractor or a long-lived transient state. 

The instability of the homogeneous state thus leads to 
bundle shortening. It is triggered by an increase of a 
and occurs for all values of ft, which suggests that the 
interaction of parallel filaments is responsible for bundle 
shortening. On the other hand, the interaction of an- 
tiparallcl filaments (ft ^ 0) induces currents which in the 
unstable regime lead to a separation of plus- and minus- 
filaments [ydj . The role of interactions between parallel 
filaments for bundle shortening can be demonstrated ex- 
plicitly for the simple case D — for which the uniform 
state is linearly unstable for all a > 0. In this case, using 
Eqs. (1) and (2), we can show that for ft = and a > 
the variance decreases monotonically, i.e., d/dt(x 2 ) < 0, 
which implies bundle shortening. 

Can the active filament interactions which lead to bun- 
dle shortening generate tensile stresses that can be used 
to perform mechanical work during shortening? In the 
bundle, tension arises due to forces generated by motors 
and due to hydro dynamic forces exerted on the filaments 
by the surrounding fluid. Let us estimate this tension a 
which has units of force. For simplicity, we take only 
local friction into account and ignore interactions be- 
tween different filaments by friction forces. Consider a 



rigid rod-like filament moving in a viscous environment 
in the direction of its axis with a friction coefficient per 
unit length rj. If this motion is induced by a locally ap- 
plied force, then the tension profile along the filament 
is piecewise linear, see Fig. 1(c). Tension is positive in 
the rear part of the filament being "pulled" and nega- 
tive in the front. Using this tension profile, we calcu- 
late an average tension profile along a filament at posi- 
tion x by summing over all contributions due to interac- 
tions with filaments lying in the interval [x — £, x + £] . 
We then calculate the tension in the bundle at posi- 
tion y, by summing the average tension of all filaments 
with centers within the interval [y — l/2,y + l/2\. Fi- 
nally, we coarse grain over one filament length and obtain 



+ + 



(y) + (i + (y)+0 + (y) + cr (y), where 



v+e/2 e 

a ±± =fj£a J dx J d^c ± (x + e)c ± (x) 

y-t/2 -I 

v+i/2 e 



T¥ft dx di [c T (x + £,)-c T (x- £)] ^(x) 



v -e/2 o 



Here, fj, f] are effective frictions per unit length which re- 
sult from the coarse-graining and differ from rj by dimen- 
sionless geometric factors. Note, that the final coarse- 
grained tension profile is independent of microscopic de- 
tails of the interaction mechanism. 

As an example, consider a homogeneous ring, i.e., peri- 
odic boundary conditions. According to the above Equa- 
tion it generates a tension a = 2afj£ 3 (cq + Cq ) , which 
for a > is positive and could then lead to active con- 
traction of the ring. This tension persists only if the 
homogeneous state is stable. As soon as an instability 
occurs, i.e., for a > a c , the ring ruptures and the remain- 
ing bundle subsequently shortens. The shortening of the 
unstable bundle generates transient inhomogeneous ten- 
sion profiles. Fig. 4 displays instantaneous density and 
tension profiles of a filament bundle during shortening. 
The density profile is symmetric, vanishing at the bundle 
ends. The tension profile has the same qualitative shape. 
The inhomogeneous tension leads to a filament current 
J(x) which is antisymmetric corresponding to filament 
motion towards the bundle center, consistent with bun- 
dle shortening. 

The effective model discussed above does not spec- 
ify the microscopic origin of active filament interactions 
characterized by a and ft. In particular, the interaction 
between parallel filaments could seem surprising jlj] . In 
order to demonstrate that such interactions can emerge 
naturally from simple motor-filament interactions we 
have performed computer simulations of a more micro- 
scopic model. We consider rigid filaments aligned 
along the x-axis. During each time-step, a motor com- 
plex creates a mobile cross-link between a randomly cho- 
sen pair of filaments. The motors are displaced towards 
the plus ends of the filaments with velocity v, possibly 
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leading to a relative filament displacement. This gen- 
erates an interaction between antiparallel filaments. A 
relative displacement between parallel filaments occurs if 
we furthermore assume that the filament ends have dif- 
ferent properties than the bulk of the filaments, namely, 
a motor which binds to or arrives at the plus end of a fil- 
ament stays attached for some time, see Fig. 1(b). This 
effect generates an effective interaction with a > be- 
tween parallel filaments ]l5| ). We have added diffusion 
steps in the simulations with D ^ 0. Our simulations 
show the same qualitative behaviors discussed above for 
the phcnomcnological model: an instability of the homo- 
geneous state which occurs below a characteristic value 
of D. In the unstable state we observe bundle short- 
ening due to the interaction of parallel filaments, inter- 
action of anti-parallel filaments leads to a separation of 
plus- and minus-filaments. Using piecewise linear tension 
profiles for moving filaments (corresponding to local fric- 
tion), we determine a(x) in our simulation by averaging 
over many numerical steps. For the homogeneous state 
we find a tension which fluctuates around a constant pos- 
itive value. The simulated density profile, tension profile 
and averaged filament current of a shortening bundle are 
displayed in Fig. 4. 

Concluding, we have demonstrated that a bundle of 
aligned filaments and motors may contract and gener- 
ate tension even if it is lacking a spatial organization. 
Using symmetry arguments, we have shown under some 
simplifying assumptions that the interaction of parallel 
filaments is important. Whether stress fibers in cells use 
the physical mechanism proposed here for contraction re- 
mains to be tested. In particular, experiments which fo- 
cus on the interaction between parallel filaments should 
be performed. Stress fibers are likely to be more complex 
than the system studied here and might also use other 
mechanisms for contraction. For example, they could 
have some yet unknown spatial organization that allows 



c) 




d) 



them to contract more efficiently, or other components 
in addition to motors and filaments might be involved. 
Furthermore, they might work close to the percolation 
threshold where large filament aggregates occur and 
which is not captured by our model. Nevertheless, we 
think that our model helps clarifying the role of motor- 
filament interactions and self-organization phenomena 
for force generation in biological cells. 
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FIG. 4. Instantaneous filament distributions (full lines, 
dots), the tension profile (dashed lines), and the current 
(dash-dotted lines) in arbitrary units of contracting filament 
bundles, (a), (b) Numerical solutions of Eqs. (1) and (2) with 
a = 0.5, ft = 0, c = 0.5, 8c = 1, and L = 10. (c), (d) Simula- 
tion of TV = 1000 filaments as described in the text, (a) and 
(c) are for earlier, (b) and (d) for later times. 
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